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Abstract 

A reduction of the Maxwell's system to a Fredholm second-kind integral equation with 
weakly singular kernel is given for electromagnetic (EM) wave scattering by one and 
many small bodies. This equation is solved asymptotically as the characteristic size 
of the bodies tends to zero. The technique developed is used for solving the many- 
body EM wave scattering problem by rigorously reducing it to solving linear algebraic 
systems, completely bypassing the usage of integral equations. An equation is derived 
for the effective field in the medium, in which many small particles are embedded. A 
method for creating a desired refraction coefficient is outlined. 

0MSC: 35J10, 70F10, 74J25, 81U40, 81V05 

PACS: 43.20.-hg, 62.40.-Hd, 78.20.-e 

1 Introduction 

There is a large literature on electromagnetic wave (EM) scattering ([!]- [3]), to name a few 
books. In t5j(see also [6], [7] and references therein), wave scattering theory is developed 
for small bodies of arbitrary shapes. The ideas and methods from [5] were used in recent 
papers [8], [11], where a many-body scattering problem was solved asymptotically for 
scalar wave scattering, as the characteristic size a of small bodies tends to zero. One of 
the applications of this theory in acoustics is a method for creating new materials with the 
desired refraction coefficients. In particular, the new materials may have a desired wave- 
focusing property: a plane wave incident at a fixed direction with a fixed wavenumber will 
be scattered by the new material with a desired radiation pattern. Another application 
(|13j. (|14j)) is a method for creating materials with negative refraction, that is, materials 
in which the group velocity is directed opposite to the phase velocity. 

In [15 the new materials are created by embedding small particles into a given material 
in a bounded domain D so that the number of the particles embedded into a small subcube 
A C -D, is given by the formula 

M{^) = ^ f N{x)dx[l + o{l)], a^O, (1) 



^key words: electromagnetic waves, wave scattering by small bodies, many-body scattering, "smart" 
materials 
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where a is the radius of the embedded spherical particles Dm = B{xm, a), Xm is the center 
of the ball B{x mjfi), o is its radius, < x < 1 is a parameter, N(^x) > is an arbitrary 
continuous in D function, which we can choose as we want, the boundary condition on 
the surface Sm of the m-th particle Dm is: 



du h{xm) ^ 

9N = ^^ on 5™, (2) 



where A'^ is the unit normal to Sm, pointing out of Dm, and h{x) is an arbitrary continuous 
function which we can choose as we want. 

The purpose of this paper is to develop a theory for the electromagnetic (EM) wave 
scattering, similar to the one developed in [8]- [12], [14], [15], in the sense that it yields a 
method for creating materials with a desired refraction coefficient. Condition ([2]) is not 
used in this paper, but formula (f73l) of this paper is a generalization of ([I]). 

In Section [2] we present a theory for EM wave scattering by a single small body. In 
Section [3] many-body scattering problem is solved. In Section [4] the negative refraction is 
discussed. In Section [5] the conclusions are formulated. 



2 EM wave scattering by a single small body 

Let D be the body and a := ^ diamZ). Assume that 

ka < 1, (3) 

where A; > is the wavenumber. The governing equations are: 

VxE = iujfiH, V X H = -iuje'{x)E inR^ (4) 

where w > is the frequency, = const is the magnetic constant, €'{x) = e > in 
D' = R^\D, e'{x) = e{x)+i^; a{x) > 0, e'{x) / OVx G R-^, e'(x) G C^{R^) is a twice 
continuously differentiable function, and cr{x) = in D' is the conductivity. From (4) one 
gets 

V xV X E = K'^(x)E, H = ^ ^ ^ , (5) 

K^{x) := LO^e'{x)fi. (6) 
We seek the solution of the equation 

V xV X E = K^{x)E (7) 

satisfying the radiation condition: 

E{x) = Eo{x) + v, (8) 

where Eq{x) is the plane wave, 

Eo{x)=£e"'^-^, k = -, (9) 
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where c = oj^jejl is the wave velocity in the homogeneous medium outside D, e = const 
is the dielectric parameter in the outside region D', fi = const, a G 5^ is the incident 
direction of the plane wave, £ ■ a = 0, £ is a constant vector, and the scattered field v 
satisfies the radiation condition 

dv ,1-. 

— ikv = o(-), r = \x\ —>■ Qo, (10) 

or ^r' 

uniformly in directions P '■= 

If E is found, then the pair {E,H}, where H = solves our scattering problem. 

Our goal is to derive a Predholm second-kind integral equation for E. 
Let us rewrite equation ([7]) as 

-AE + V{V-E)-k'^E-p{x)E = 0, (11) 

where 

p{x) := K'^ix) - k^, p{x) = OmD'. (12) 
Note that ImX^(x) = Imp(x) > 0. This will be used in the proof of Claim 3 below (see 



It follows from d?]) that 

= V • {K^{x)E) = VK'^ix) • E + K^{x)V • E. (13) 
From (fT2|) and ([13]) one gets 

- - k'^E - p{x)E - V{q{x) ■E)=0, (14) 

where 

q{x):= 9(^)=0iii^- (15) 

From and ([SD one gets 

E = Eo + g{x, y) {p(,y)E{y) + V^(g(y) • E{y))^ dy, (16) 



where 



^ik\x-y\ 



9ix,y):=— :. (17) 

47r|x — y\ 

If q vanishes at the boundary of D, i.e., V-fC^(x) = Vp{x) vanishes at the boundary of D, 
then, after an integration by parts in the last term of (jl6p . and taking into account that 
Vxg{x,y) = —Vyg{x,y), one gets an equivalent equation, which we also refer to as (fT6]l : 

E = Eo+ g{x,y)p{y)E{y)dy + g{x,y)q{y) ■ E{y)dy. 

J D JD 



Since we have assumed K'^{x) ^ and K'^{x) G C^(IR'^), it follows that equation (jl6p 
is a Fredholm equation of the second kind, because the integral operator T in ()16p is 
compact in H^{D), where H^{D) is the usual Sobolev space. 
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Indeed, the operator Bf := Jj^) g{x,y)f{y)dy acts from L^{D) into H^{D), and the 
operator BiE := p{y)E + V(g(y) • acts from H^{D) into L^{D). Thus, T acts from 
H^{D) into H^{D) and is, therefore, compact by the embedding theorem. 

Let us summarize the result. 



Lemma 1 // K'^{x) G C^(M^), K'^{x) ^ 0, K'^{x) = k"^ > in the exterior domain 
D' := R'^ \ D, then the operator T in (|16p is compact in H^{D), so that equation (|16p zs 
of Fredholm type in H^{D). 

Lemma 2 Equation (jl6p is uniquely solvable in H^{D). 

Proof. Let us prove Lemma 2. It is sufficient to prove that the homogeneous version 
of equation (jl6p has only the trivial solution. If E solves the homogeneous equation (|16p . 
then E solves equation ([7]) and satisfies the radiation condition (jlOp . The only solution to 
([7]) satisfying (jlOp is the trivial solution E = 0. 
Let us give details and prove the above claims. 

Claim 1 . A solution to (fT6|) satisfies equation d?]), and, consequently, p^ . that is, 
V • = 0. 

Thus, Claim 1 states that (fTBI) is equivalent to the original equation ([7]), which is not 
at all obvious. 

Claim 2 . The only solution to (HUD is ^ = 0. 

Proof of Claim 1. HE solves (1161) and is of the form (IHl), then 



(-A - k'^)E = p{x)E + V(g(x) • E). 
Rewrite this equation as 

V X V X £; - V(V • E) - V{q{x) ■ E) = K'^{x)E, 

or 

_ _ ^ ^fK^(x)V-E + VK^(x)-E\ 

VxVxE-Vi^-±^ J =^'(^)^ (18) 

Denote 

■ {K\x)E) :=V(x). 
Taking divergence of equation (fTHj) one gets 

- AV'-i^^(x)^/; = in (19) 

The function ip satisfies the radiation condition (llOp . 

Claim 3 . The only solution to (jl9p . which satisfies condition (jlOp . is V' = 0. 

We prove Claim 3 below. Assuming that this claim is proved, we infer that = 0, so 
V • {K^{x)E) = 0, and equation holds. 

Claim 1 is proved. □ 

Let us prove Claim 3. Equation (jlOp can be written as: 

- All; - k'^ip - p{x)'il; = inM^ (20) 
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where p{x) is defined in equation (I12|) . and p{x) = in D'. It is known that if > 

is a bounded domain, and satisfies the radiation 
condition (fTOl) and equation (f20]) . then ^ = (see, e.g., [SJ). 

For convenience of the reader, we sketch the proof. From (j20p and its complex conju- 
gate one derives the relation 

ipA^ --ip Alp -2iIm.p{x)\Tp\'^ = 0. (21) 

Integrate (I2ip over a ball Bji centered at the origin of radius R, and use the Green's 
formula to get 

-tp^)ds -2i f Imp(x)\M'^dx = 0. (22) 
\x\=R dr dr' Jd 

By the radiation condition (jlOp for ip one can rewrite ()22p as 



2ifc lim 



\ip\^ds-2i Imp(x)|V'pdx = 0. (23) 



|x|=_R 

Since lmp{x) > by our assumptions, and /c > 0, it follows that 



lim / |V'pdx = 0. (24) 

R^°°J\x\=R 

Relation (j24p and the equation 

(A + k'^)'ip = |x| > fio, (25) 

where D D, imply ip = for |x| > i?o- (See |7, p. 25]). 

By the unique continuation principle for the solutions of the homogeneous Schrodinger 
equation (j20p . it folows that ip = 0. Claim 3 is proved. □ 

Let us prove Claim 2. This will complete the proof of Lemma [2j 
If E solves ([7]) and satisfies pO|) . the the pair {E,H}, where H is defined by the 
second formula in ([5]), solves the homogeneous Maxwell system satisfies the radiation 
condition, and {E,H} is in i?^^^(]R^). It is known (see, e.g., [4J) that this implies E = 
H = 0. Lemma [2] is proved. □ 

From formula (jl6p assuming that a := ^diaml? is small, /ca <C 1 and the origin is 
inside D, one gets: 



pik\x\ 1 

E = Eo + ^^ 

\x\ 47r 



as |x| — > oo, 1^ := /5, and 



r e-'>^^y{p{y)E{y) + V{q{y) • E))dy[l + 0(^)] (26) 



Thus, 



e-^^P-y = l + 0{ka), \y\<a. (27) 
E{x) = Eo{x) + / {piy)E{y) + V{q{y) • E))dy[l + 0(^)]. (28) 
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Here we have used the estimate 

(IxP - 23; • y + |yP)2 = |x| - /3 • y + O(^), |x|>a>|y|, ():=xl\x\ 

\x\ 

which yields 



e + 0(--) + 0(--- 



47r|x-y| 47r|x|(l + 0(^)) 47r 

where |y| < a < If ka < 1, then O(^) < 0(|||), and we get estimate Usuahy 
the error term in (j28p is written as 0{-^) as |x| oo, but for small a the term is 
small already for \x\ ^ a, and x may be still small. 

If D does not contain the origin and Xm G then one can rewrite ()16p as: 

(){x, Xfn 

) / Piy)E{y)dy + Vr,g{x,Xm) / q{y) ■ E{y)dy 



+ / [g{x,y) - 9{x,Xm)\p{y)E{y)dy + Vx I [gix,y) - g{x, Xm)]qiy) ■ E{y)dy, 
Jd Jd 

(29) 

where 5 is defined in (jl7p . Let 

K^ := / piy)E{y)dy, f q{y) ■ E{y)dy, (30) 

Jd Jd 

where p{x) is defined in (jl2p and g(x) is defined in (jl5p . 
If ^ 

d := |x — Xm| ^ a, a := - diamZ), 

then 

/t 1 

15(2;,^) - 9{x,Xm)\ = \Vyg{x,y) ■ {y - Xm)\ < camax(-,^), (31) 

where c > stand for various constants independent of a, y is an intermediate point, and 
d = \x — Xm\- Thus, 



k 1 

< ca max{-, -^)\p{xm)E{xm)\\D\, 



[g{x,y) - g{x,Xrr,)]p{y)E{y)dy 

ID 

where \D\ = 0{a'^) is the volume of D, and we assume that p{x) and E{x) are continuous 
in D. 

On the other hand. 



g{x, XfYi ) / p{v)E{y)dy 



D 



^(Mi^I!i)MMM)>>j„, (32) 



provided that § ^ 1 and ka <^ 1. These assumptions hold throughout the paper. 
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We have checked that 
Similarly, 
where 



J'm ^ I'm^ a d, ka '^1, 



r :-- 



Vr,g{x,Xm) q{y) ■ E{y)dy 



D 



J' - 



Va. / [g{x,y) - g{x,Xra)]q{y) ■ E{y)dy 



Also, 



D 



(33) 

(34) 
(35) 



Therefore, in the region \x — Xm\ ^ ^ one can calculate E{x) by the approximate 
formula 

(36) 

with the error 0(^ + ka). Here Vm and i^m are defined in ([30|) . 
Let us derive formulas for the vector Vm and the scalar Vm- 
Multiply by p{x) and integrate over D to get: 



Vjn = Voni+ j p{x)g{x,Xm)dxVm+ / p{x)V xg{x , Xm)dx 
JD Jd 



(37) 



where 



Vom ■■= / p{x)Eo{x)dx, 



D 



and we have neglected the two last terms in (j29p . These terms are smaller that the ones 
we kept, as we prove soon. 
Let 

p{x)g{x,Xm)dx := am, / p{x)Vxg{x,Xm)dx := Am- (38) 



D 



D 



Take dot product of (j29p with q{x), integrate over D and again neglect two last terms to 
get 

(39) 



where 



^Om ■= / q{x) ■ Eo{x)dx, 
Jd 



Br, 



:= / q{x)g{x,Xm)d 
Jd 



D 

X, bm 



D 



l{x) ■ Vxg{x,Xm)dx. 



From (|37j) one gets 



Vm. = h V„ 



1 O^m. 1 ^r. 



(40) 
(41) 

(42) 
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Note that if the body D is small, then \am\ < 1, and actually \am\ ^ 1, so (I42p is well 
defined. Also, |6m| < 1 if is small. From (|42p and (j39p one obtains: 

_ (1 — am)t'Om + • Vbm ^^^^ 



(1 (2^1 )( 1 &m ) Bqji ■ AfYi 

and 

Y _ ^Om _|_ (1 ~ O'mj^Om + -^m " ^Om (44) 

1 '^m 1 (1 <^m)(l ^m) -f^m ' 

Formulas ([36]) . ([l3j) . and (jUj) sofce t/ie scattering problem of EM waves by one small body 
D located so that Xm G D. 

Let us now check that one may neglect the terms we have neglected. 

For example, one has 

dxp{x) / [g{x,y) - g{x,Xm)]p{y)E{y)dy = / dyp{y)E{y) / dxp{x)[g{x,y) g{x, Xjji 
Jd Jd Jd 

dyp{y)E{y) [V'i(y) - ^i(xm)] 
< / dy\p{y)E{y% 



D 



D 



D 

because that V'i(y) '■= fodxp(x)9(X'y) ^ continuous function. 

The function 'ijji{y) is continuous if p{x) G L?'{D), since ipi G i/;^Q^(M^) and the embed- 
ding from H^{D) into C{D) is continuous (and compact if D C is a bounded domain, 
which we have assumed). 

A similar argument holds for ■(/'2(y) := [^x9{x,y)— Vx{9{x,Xm.)]-'l{y)dy Si,nd iustiG.es 
the possibility to neglect the second term we have neglected. 

To summarize: 

If p(x) and q{x) are continuous functions then the scattering of EM waves by a small 
body D is solved analytically by the formulas (j36p . (j43p . (|44p . with an error of order 

3 EM wave scattering by many small bodies. Creating ma- 
terials with a desired radiation coefficient. 

The basic equation now has the form 

M ^ M 



E{x) = Eo{x) + g{x,y)p{y)E{y)dy +^Vx j g{x,y)q{y) ■ E{y)dy 

M M 
= Eo{x) + ^ [g{x,Xm)Vm + Vxg{x,Xm)l^m\ + 



m=l m=l 

where M is the number of small bodies. 



Jm-= I [gix,y) - g{x,Xra)]piy)E{y)dy, (46) 



D„ 
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and 

Km:=V^ [g{x,y)-g{x,x^)]q{y)-E{y)dy. (47) 

J Dm 

We assume throughout, that the points £ are distributed in D so that their 
number AA(A) in any subdomain A E D of is given by formula (j73p (see below) with 
'^{a) = c?~^ ^ X > (see Lemma 3 below). As a ^ and M oo, the set of the points 
Xm becomes everywhere dense in D. 

Neglecting Jm and Km,-, let us derive a linear algebraic system (las) for finding Vm and 
Vra- If Kra and fm, 1 < ^ are found, then the EM wave scattering problem for M 
small bodies is solved by the formula 

M 

E{x) = Eo{x) + [ 

(48) 

with an error 0(^ + ka) in the region mini<m<A/ 1^; — Xm\ := d ^ a. 

To derive a linear algebraic system for Vm and Um, multiply (|15]) by p{x), integrate 
over Dj, and neglect the terms Jm and Km, to get: 

M 

= % + Yl i^jrnVm + Bjml^m) , 1 < J < M, (49) 
m=l 

where 



Vbj := / pix)EQ{x)dx, ajm= / p{x)g{x, Xm)dx, Bjm := / p{x)Vxg{x,Xm)dx. 

JDj J Dj J Dj 

(50) 

Take the dot product of (PS]) with g(x), integrate over Dj, and neglect Km, Jm, to get: 

l^j = VQj + ^ [CjmVm + djml^m) , (51) 
m=l 

where 

z^Oj = / q{x) ■ Eo{x)dx, Cjm-= / q{x)g{x,Xm)dx, djm ■= / q{x) -V xg{x,Xm)dx. 
JDj JDj JDj 

(52) 

Equations (j49p and (j5ip form a linear algebraic system for finding Vm and , 1 ^ ^ . 
This linear algebraic system is uniquely solvable if A;a ^ 1 and a <^ d. The elements Cjm 
and Bjm are vectors, and ajm, djm are scalars. Under the conditions 

M 

~ ~ m=l 

one can solve linear algebraic system ()49p and (j5ip by iterations. In (j53p ||-Bjm|| and ||Cjm|| 
are the lengths of the corresponding vectors. Condition (|53p holds if a ^ 1 and M is not 
growing too fast as a ^ 0, not faster than 0(a^'^). 
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Consider the limit, as a ^ and M — > oo, of the field E{x) := E[x, a). Let us assume 
that the number AA(A) of small particles in any open subset /\ <Z D \s given by the 
formula: 

AA(A) = [ N{x)dx [1 + o(l)] , a ^ (54) 



where < xi < 1 is some number, and N{x) > is a continuous in D function. If ([54 
holds, then the distance between two neighboring small particles is of the order 0{a^^). 
We choose xi = (3 — x)/3 above formula ([75]) . 
Let us write (HSl) as 



E{x) = Eq{x) + g{x,Xm) p{y)dyE{xm) + "^x9{x,Xm) I \ " ^iy)'^y 

(55) 

where the terms Jm and Km are neglected. 

Choose the function p so that suppp = Um=i ^m, where Dm := Bm '■= {y '■ \y — Xm\ < 
a}, that is, the small particles are balls of radius a centered at the points Xm, and suppp 
is the support of p. Choose p of the form: 

p{y)=p{r,a), y £ Bm, r := \y - Xm\, I < m < M, (56) 

where the function p{r, a) should be chosen so that 

/ p{y)dy = 0{a^), a ^ 0, s := 3xi = 3 — x, (57) 
where < x < 3. One may like to have Zm = 0{a^) as a ^ 0, where 

with the same s as in (j58p . If this would be possible, then it would make possible passing 
to the limit as a — > similarly to the passing done in [T3]. However, as we will see below, 
the relations (j57p and Zm = 0{a^) are not compatible, in general. 
One has: 

V7^\ '/ \ y ~ ■^'m ! dp , , 

Vp{y) = p {r,a)r", := , p ■= —. (59) 

r ar 

Here is a unit vector and r is the length of the vector y — Xm- Let 

T?t \ TPl \ I dEj{Xm) I ,0/' 2\ Vi ~ {Xm)i /„„n 

E{y) = E{xm) + rcj — ^ + 0{a ), r"^ := , (60) 



{xm)i ■= Xm • Cj, a • 6 is the dot product of vectors a and 6, ei ■ Cj = 5ij is the standard 

0, i + j, 



Cartesian basis of M'^, and (5; 



Over the repeated indices in ([60]) and everywhere below summation is understood. 
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If p{y) = p{r, a) then 



m ) 

Bri 



A;2 +p{y) 



dy = (61) 



for any constant vector E{xm) by symmetry. 
Thus 



By symmetry the integral in (1621) vanishes if i ^ j. This integral for z = j is equal to 
drr^rp'{r,a) r '"^'^ 



^i-= I To / d^sin^ / dv^r^, i = l,2,3, := r^. (63) 

'o +p{r,a) Jo 







One has r| = cos^ 0, = sin^ sin^ ip, rf = sin^ cos^ ip, so 



27r 4^ 

dO sin 9 j diprf = — , 
3 



and, with r = ta, one gets: 



Y^ = ^-la^ r ffy;!'i +o(a^), a^O, ^ = 1,2,3. (64) 

Thus 

Z™~^a4/(a)V-i?(x™), lia) := [' r = ta, (65) 
3 Jo k^+p{at,a) 



and 



jm '■= / pdy = 4:t:o? I dtt^p[at,a) := Ana^ j{a). 
JBm Jo 



Let 



7m 



(1 - t)2/i(t), < t < 1, t := ^, X = const > 0, 



(66) 



p(r)=p(r,a) = | v j k " ^ «' ^ ' (67) 

where 7m, Im7m > 0, is a number we can choose at will. Then 

f /"^ 

jm ■■= / Piy)dy = — Tm / (1 - tfh{t)t^dt := cvma^-'', (68) 
JB„ a"" Jo 



where cim is proportional to 7^, so we can also choose cim at will. By (j65p . 

Jo I^ + 7m(l-t)2/l(t) ^ ^ 

Can one choose h{t) so that /(a) = Co^^[l + o(l)] as a ^ 0? If yes, then Zm and are 
of the same order as a ^ 0. 
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In fact, it is not possible to choose such an h{t). One has 

I{a) = t^ ln[zy + 7^(1 - tfhit)] In - 3 / t^ ln[u + 7m(l - tfh{t)]dt. (70) 

Jo 

Thus 

1(a) = In i/ + 0(l) = xln a + 0(1), a ^ 0. (71) 

Therefore no choice of /i G C[0, 1] can ensure the asymptotics I{a) = Ca~^[l + o(l)] 
as a — > 0. 

We have obtained the asymptotic formulas: 

jm = ci„a3-^ Zm = ■ E)ixm)a^ In a[l + o{l)], a ^ 0. (72) 

Let us take h{t) = 1. Then cim = 7m/30. To pass to the hmit a ^ in formula ([55]) we 
use the following lemma, which a generalization of a lemma from [8] : 

Lemma 3 // / € C(Z)) and distributed so that 

MiA) = ^- [ N{x)dx\l + o{l)], a^O, (73) 
fW J A 

for any subdomain A C D, where (p{a) > is a continuous, monotone, strictly growing 
function, ip{0) = 0, then 

lim Tf{xmMa)= [ f{x)N{x)dx. (74) 
-^^^ Jd 

Proof. Let D = UpAp be a partition of D into a union of small cubes Ap, having no 
common interior points. Let |Ap| denote the volume of Ap, 6 := maxpdiam Ap, and y be 
the center of the cube Ap. One has 



limV/(y(P))iV(y(P))|Ap| = / f{x)Nix)dx. 



The last equality holds since the preceding sum is a Riemannian sum for the continuous 
function f{x)N{x) in the bounded domain D. Thus, Lemma 3 is proved. □ 

Let (p{a) = a^-^, x > 0. Then, by formulas ([58]), dSSD and ([72]), one gets: 

hm -^P- = 0, lim = cim / 0. (75) 
a^O (p[a) a^O (p[a) 

Therefore, by Lemma El the right side of formula (1550 has a limit as a — > provided 
that (j54p holds with 3xi = 3 — x, where x is defined in formula (j68p . The limiting equation 
(I55p takes the form 

Ee{x) = Eo{x) + [ g{x,y)C{y)Ee{y)dy, (76) 
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where Ee(x) := lirria^o E{x) is the effective field in the medium when the number of small 
particles tends to infinity while their size a tends to zero, and 

C{Xm) ■■= CimN{Xm). (77) 

Formula (j77p defines uniquely a continuous function C{x) since the points 
distributed everywhere dense in as a — > 0. The functions C(x) can be created at our 
will, since it is defined by the numbers cim and N{xm), which are at our disposal. 

Apply the operator + k'^ to (f76]) and get 

[V^ + IC\x)]Ee = 0, /C2(x) := A:2 + C7(x) :=yfcV(x). (78) 

Therefore the (refraction) coefficient n'^(x) is defined by the formula: 

n\x) = l + k~^C{x). (79) 

The function IC'^{x) is quite different from the function i^^(x), defined in ([6j). The function 
C(x) and, therefore, n^(x) depend on the choice of A^(x) in (j73p and on the choice of the 
coefficients cim in ([72]) . The cim is determined by 7m, defined in ([67|). Since one is free to 
choose 7m and N{x) as one wishes, one can create a desired function n^(x) by choosing 
suitable 7m and N{x). 

Note that the total volume of the small particles is negligible as a — > 0. Indeed, this 
volume equals to the volume of one particle, which is times the total number Af{D) 
of the smah particles in D. One has x > 0, so Af{D)^a^ = 0{a-^^-'^^a^) = 0{a'^) 
as a ^ 0. 

The limiting effective field Eg satisfies equations ([76|) and ([78]) . but does not, in general, 
satisfy the equation 

V ■Ee = 0. (80) 
Indeed, the solution to ()76p is unique, and it follows from ()78|) that 

V2r/ + /C2r/ = -V/C2(x) -Se, iC^ (x) := k"^ + C (x) , (81) 

where rj := V ■ E^. Since (|8ip is inhomogeneous, and V/C^ • E^ ^ 0, in general, its solution 
r/ 7^ 0, in general. 

The unique solution to ([76]) does not satisfy, in general, the Maxwell's equation ([7]). 
Thus, in the limit a — > one gets a vector field E^ which satisfies a new equation (|76p . 

What is the physical meaning of this new equation ()76p ? This equation is equivalent 
to equation (j78p the solution of which is required to satisfy the radiation condition. If one 
rewrites equation ([78ll in the form: 

V XV X Ee= IC'^{x)Ee + VV ■ Ee, (82) 

then one can interpret the term VV • Ee as the "current" itofiJ in the Maxwell's equation. 
Indeed, Maxwell's equations 

V X E = iuj^iH, V X H = -iuje'{x)E + J, 
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imply 

V xV X E = K'^{x)E + iujfiJ. (83) 

Therefore the term {iujfi)^^W • := J can be interpreted as a "non-classical" current: 
the classical current is of the form aE, where a = aij is the conductivity tensor. In our 
interpretation tensor cr is a differential operator: a = Cij = q^.q^ . , so aE = VV • E. 

A different interpretation of this term may be based on the assumption that the vector 
of the induction De in the limiting medium is related to Eg non-locally: 

Deix) = eEeix) + / x{x,y)Ee{y)dy, 
Jd 

where xi^^u)-, the susceptibility tensorial kernel, is the kernel of the non-local operator 
defining the polarization vector in the limiting medium. Such non-local operator of po- 
larization will produce the term VV • Ef, if x(x, y) is a distributional kernel of the form 
X(x,y) = (u;2pt)-iVa.(5(a; - y)Vy. 

Indeed, Maxwell's equations in this case are: 

V X E = iujfj.H, V X H = —iioeE — iuj j x{x-, y)E{y)dy. 

JD 

With the above choice of x{Xjy)i these equations imply: 

V xV xE = IC'^{x)E + VV-E. (84) 

This is equation (|82p. which is equivalent to (|78p . The function e in the above Maxwell's 
equations is related to /C^ by the usual formula io'^fie = K? . 

It would be of interest to understand the properties of the limiting medium from an 
experimental point of view. 



4 Spatial dispersion and negative refraction 

If v?{x) = n?{x,uj), so that uj = ci;(K), then one says that spatial dispersion takes place. 
Here K is the wave vector, |K| = K. Its direction is given by the unit vector K*^, which 
is oriented along the direction of the propagation of the wave front. The group velocity 
Vg is defined by the formula 

Vg = Vkw(K), 
while the phase velocity Vph is defined by the formula 

^ ^^K", K0:=-, i^:=|K|. 

Vh |j^| ' ^' II 

Negative refraction means that the group velocity in a medium is directed opposite to the 
phase velocity. 

A sufficient condition for an isotropic material to have negative refraction was derived 
in [13]. Namely, in an isotropic medium uj = a;(|K|), and 

ujn{x,uj) = cK, K = |K|, 
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where c > is the wave speed in the free space. Applying the gradient with respect to K 
to the above equation yields 



If lmn{x,uj) = 0, then (j85p shows that V-kujCK) is directed opposite to K*^ if and only if 



To create a material with negative refraction it is sufficient to create an isotropic ma- 
terial with a refraction coefficient n{x,uj) satisfying inequality (j86p . This can be achieved 
by choosing the constants cim in (f68]) depending on uj so that the function C{x,uj) would 
generate by the formula 



the refraction coefficient which is a real- valued function satisfying inequality (j87p . 

5 Conclusions 

A theory of electromagnetic (EM) wave scattering by many small bodies is developed. 
A rigorous reduction of the many-body scattering problem to solving linear algebraic 
system is given bypassing numerical solution of the integral equations. An equation for 
the effective field in a medium consisting of many small particles is derived. It is shown 
that by choosing large number of small inhomogeneities with special properties one can 
create a medium with a desired refraction coefficient. 
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